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We present a framework to elucidate the existence of accidental contacts of energy bands, par- 
ticularly those called Dirac points which are the point contacts with linear energy dispersions in 
their vicinity. A generalized von-Neumann-Wigner theorem we propose here gives the number of 
constraints on the lattice necessary to have contacts without fine tuning of lattice parameters. By 
counting this number, one could quest for the candidate of Dirac systems without solving the sec- 
ular equation. The constraints can be provided by any kinds of symmetry present in the system. 
The theory also enables the analytical determination of fc-point having accidental contact by selec- 
tively picking up only the degenerate solution of the secular equation. By using these frameworks, 
we demonstrate that the Dirac points are feasible in various two-dimensional lattices, e.g. the 
anisotropic Kagome lattice under inversion symmetry is found to have contacts over the whole lat- 
tice parameter space. Spin-dependent cases, such as the spin-density-wave state in LaOFeAs with 
reflection symmetry, are also dealt with in the present scheme. 
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I. INTRODUCTION 

An issue regarding the contacts of energy bands has 
long been studied from the early stage of sohd state 
physics^. The one currently attracting wide interest 
is called Dirac point which is characterized by the lin- 
ear splitting of energy bands in its vicinity. This dis- 
tinguishing band structure leads to the exotic phenom- 
ena in the electronic transport, orbital diamagnetism, 
etc^. The Dirac points have been observed and stud- 
ied in numbers of materials: graphene^^^, organic solids 
such as a-ET2l3^'^, LaOFeAs*'^ in its spin-density-wave 
state, and those with spin-orbit interactions^''. The 
ones in graphene are located at K and K'-points at 
the corners of the hexagonal Brillouin zone^^'^^. How- 
ever, in a-ET2l3, they fall on some general fc-points 
inside the Brillouin zone^'^'^^. In two-dimensional(2D) 
models, their presence are shown in several examples, 
such as the isotropic^ ^ and anisotropic^^ honeycomb 
lattices, the anisotropic square lattice under the pres- 
ence of a TT-flux^^, the non-Bravais anisotropic square 
ones^^, and the isotropic Kagome lattice^'', etc. The 
first-principle^^, semi-empirical calculations^^, and per- 
turbative calculations^^ are used to show the existence 
of Dirac points in Q;-ET2l3. However, almost all these 
previous studies focus only on the specific and detailed 
models or materials, and lack the general viewpoint of 
how to design the systems with Dirac points in their en- 
ergy hands (Dirac systems). 

The contact, or equivalently, degeneracy of energy 
bands at the same fc-point, is classified into the essential 
or accidental (nonessential) ones, according to whether 
or not we can specify that fc-point in advance without 
solving secular equation. By definition, the essential de- 
generacy is forced to take place at special k-points, while 
the accidental ones can sometimes occur at general k- 
points. Thus, the Dirac points found in Q;-ET2l3 is clearly 



accidental. Although the general theory on the essential 
degeneracy is well established, that on accidental one has 
been only poorly explored. Albeit, the demand for the 
latter theory is developed recently, with a need to design 
Dirac systems. 

In the present paper, we propose a feasibility (gener- 
alized von-Neumann) theorem, which provides a general 
treatment on the accidental degeneracy of energy bands. 
On the basis of this theorem, one could explore the can- 
didates of Dirac systems only by counting the number 
of constraints on the lattice, without solving the secular 
equation. Our framework also provides a practical pro- 
cedure to search for the fc-points at which the accidental 
degeneracy takes place (see the introductory part of §111). 
By this procedure, one can selectively pick up only the 
degenerate solutions of the secular equation, which is in- 
dispensable in the analysis of the accidental degeneracy. 

In §11, we present a feasibility theorem in detail. In 
§11 F, the relation between the symmetry of the sys- 
tem and the number of constraints is discussed. Then, 
we demonstrate in §111 that Dirac systems are easily 
designed on various 2D lattices. In fact, the present 
study lead to the rediscovery of another Dirac sys- 
tems in multi-band two-dimensional organic crystal, 
(DIED0)2X (A=Cl,Br)20. Those who are not interested 
in the details of the formulation can first go through §11 A 
and then directly apply the practical procedure explained 
in the first part of §111. 



II. GENERAL CONSIDERATION 
A. Preliminary definitions 

We consider a lattice with orbitals per primitive unit 
cell. Our primary concern is the non-interacting system, 
while interaction effects can be treated via the one-body 



2 



approximation. The effect of external magnetic field is 
also not considered, so that the system is invariant under 
the time-reversal operation. However, our formalism can 
be easily extended to the systems in a uniform magnetic 
field, by the introduction of an extended magnetic unit 
cell with an integer magnetic flux quanta and a corre- 
sponding folded magnetic Brillouin zone. Here, the spin 
degrees of freedom is also formally neglected, whereas we 
can take it into account by regarding up- and down-spins 
as different orbitals as mentioned in §H G, which simply 
doubles the number of the orbitals. 

Following Ref. 1, we are using the term contact to de- 
note the degeneracy of energy bands at the same fc-point. 
Let the contact be called rid-dimensional, when it takes 
place on a nd-dimensional manifold in the fe-space. In 
particular, the zero- and onc-dimcnsional ones arc named 
point contact and line contact, respectively. The point 
contact at fc = fco is designated as Dirac point, when two 
bands split linearly in energy for any small Sk = k — kg 
and form an elliptic Dirac cone. Besides, we can also 
consider the situation for rig > 2, where three or more 
bands are touched together at the same fc-point. It is 
expressed as m-fold contact, where m is the number of 
bands touched. 

The lattice system is characterized by a set of param- 
eters, which we call lattice parameters, t = (ti,i27''')- 
For example, in the tight binding models, t indicates the 
transfer integrals, the energy levels of the sites, and the 
spatial coordinates of sites, and so on. In the first prin- 
ciple calculations, t denotes the Fourier components of 
the self-consistent lattice potential. If the contact takes 
place in finite and connected region in the lattice param- 
eter space, the contacts are called feasible. By contrast, if 
its presence needs fine tuning of the lattice parameters'^, 
or equivalently, if it occurs only at an isolated point in the 
lattice parameter space, the contact is considered to be 
unfeasible. The latter is unworthy of attention, because 
it is too fragile to be observed in the realistic materials. 

As will be eventually shown in §11 C, in order to ex- 
plore the candidates of lattices with feasible contacts, one 
has only to count the number of constraints on the lat- 
tice parameters, Uc, which are often provided by symme- 
tries. Actually, the dimension, n^, of the feasible m-fold 
contact should satisfy the feasibility ( or generalized von- 
Neumann-Wigner) condition, 

n-d = — rn^ -|- 1 + ric > 0, (1) 

where denote the number of unknown variables used in 
searching for the contact. The number of constraints and 
unknowns satisfy < Uc < — 1, and < < d, re- 
spectively, where d denotes the spatial dimension. When 
one searches general k-points on the Brilloiun zone, the 
number of unknowns is given as riu = d. Sometimes, the 
number of constraint, Uc, increases at special k-points. 
In this case, the number of unknowns, riu, should be re- 
defined as the dimension of manifold of these special fc- 
points. As an example, let us consider a two-dimensional 
lattice {d = 2) with a reflection symmetry. The fc-points 



on the symmetry axis must be distinguished from the 
general fe-points, since one finds larger number of con- 
straints. The number of unknowns is reduced to = 1 
there. 

In our context, the contact is considered to be essen- 
tial, when there exist the maximum number [n^ = in^ — l) 
of constraints on some special fc-points. Actually, one al- 
ways finds nu-dimensional feasible contact on these spe- 
cial fc-points. Otherwise, the contact is accidental, and 
thus we need to search over the fc-space under the prop- 
erly chosen lattice parameters. It is noteworthy that the 
accidental contact can be feasible on the general fc-points 
in contrast to the essential one. 

Among the feasible contacts, the one at fc = fco is con- 
sidered to be stable, if we can find a new contact in the 
neighborhood of fe = fco after any infinitesimally small 
variation of lattice parameters, t — > t + 5t under the 
above mentioned constraints. Conversely, contacts be- 
come unstable, just when they are created or destructed 
by some variation of lattice parameters. More intuitively, 
the stable contacts are realized inside the lattice param- 
eter region showing contacts, while the unstable ones on 
its edge. 

B. Contacts in two- band systems in 2D 

Let us begin with the two-band systems (ng — 2) in two 
dimension. The properties of point contacts in these sys- 
tems have already been investigated in detail, and most 
of such studies depend on the topological arguments^'''. 
In this section, we present some elementary explanation 
as an introduction to the multi-band cases. 

The Hamiltonian at the Bloch wave vector, fc = 
{ki,k2), is expressed by a 2 x 2 Hermite matrix. Fol- 
lowing the notations of Ref. 22, we expand it as 

H{k) = Eo{k)i + R{k)-a, (2) 

where i?o(fc) and R{k) — (i?i(fc), i?2(fc), i?3(fc)) are real 
functions of fc, and a- = (tJi , 02 , (Tg) are the Pauli matrices 
defined by 

<^i= (? J), <^2= (° ^^^{l \)- (3) 

Such an expansion is always possible, because the unit 
matrix, /, and the Pauli matrices form a complete basis 
for the 2x2 Hermite matrices. Regarding the second 
term in Eq. (2) as a fictitious Zceman splitting, we im- 
mediately obtain the eigcnencrgies as 

e±(fc) = i?o(fc)±|i2(fc)|, (4) 

which shows that a contact takes place at fc = fco, where 

R{k) ^ (5) 

is satisfied, and its energy is given as eg = Eo^ko). It is 
intrinsic that three conditions (i?i = i?2 = ^3 = 0) must 
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be simultaneously fulfilled to have a degeneracy, which is 
known as von-Neumann-Wigner theorem^'^'^*. 

Here, we consider only the point contact located at a 
general fe-point. Then, we have two unknowns, fci and fc2, 
and thus Eq. (5) is overdetermined. In order to make the 
point contact at fc = fco feasible, a constraint is required, 
which reduce the number of conditions by one. In the 
vicinity of fc = fcoj this constraint should be expressed as 



s ■ R{k) = 0, 



(6) 



within the linear order of R{k), where s is a nonzero 
three-dimensional (3D) vector. As will be mentioned in 
§11 F (in the context of Eq.(24) for multi-band systems), 
this constraint is usually attributable to the symmetries 
of the system. 

If a point contact arise at fc = fco, the energy bands 
are expanded within the linear order of (5fc = fc — fco as 



e±(fco + Sk) = eo + A- dk± \XSki + Ydk2\ 



(7) 



with 



(8) 

This contact is identified as a Dirac point, when e± form 
an elliptic Dirac cone in the vicinity of fc = fco^^i or 
equivalently, if two vectors, X and Y , are linearly inde- 
pendent. The Dirac cone tilts, ii A 0^'^. 

A feasible point contact, which is realized at k = kg 
for the lattice parameter, t = to, is classified to stable 
and unstable ones, whether or not we can find a new 
point contact in the vicinity of the original one after the 
infinitesimally small change of lattice parameter, t = to + 
6t, under the constraint of Eq. (6). Within the linear 
order of Sk and St, we obtain 



R{ko + Sk,ta + St) = R{ka - 
= XSki 



Sk,to + St) - R{kn,tn) 
-Y6k2+SRst, (9) 



where two vectors, X and Y, defined in Eq. (8) are eval- 
uated at t = to, and another vector is introduced as 



5Rst = (St ■ Vt)R\k=ko,t=to- 



(10) 



Due to the constraint of Eq. (6), three vectors, X, Y, and 
Rst are coplanar, satisfying s ■ X = s -Y = s ■ SRgt = 0. 
If two vectors, X and Y are linearly independent, the 
feasible point contact is stable. Actually, the equation, 
i?(fco + Sk, to + St) = 0, is solvable as a linear equation 
of Sk, and a new point contact is found. In other words, 
Dirac points are, by definition, stable. 

Conversely, whenever the feasible contact is unstable, 
or equivalently, just created by some lattice parameter 
change, X and Y are linearly dependent, and thus there 
exists a vector, Sq ^ 0, satisfying XSqi + YSq2 = 0. 
Thus, the band splitting A{Sk) = e+(fco + Sk) — €-{ko + 
Sk) = 2|i2(fco + Sk)\ show a quadratic dependence on Sk 
when Sk is parallel to Sq. In other words, the unstable 
point contact is given as a doubly degenerate solution of 
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FIG. 1: (Color online) (a) Schematic illustration of the merg- 
ing of two Dirac points along the direction parallel and per- 
pendicular to 5q (direction along which the Dirac points 
merge), (b) Merging of four Dirac points consisting of upper- 
and lower-pairs along the direction parallel to Sq's of each 
pair, (c) Sketches of energy bands when the band overlap is 
present (left panel) and absent (right panel). 



Eq. (5), i.e., as the consequence of the merging of two 
Dirac points, as depicted schematically in Fig. 1(a). The 
pair creation and destruction of Dirac points implies that 
there are always even number of Dirac points. This is 
nothing but Nielssen-Ninomiya (Fermion doubling) the- 
orem in 2D lattices^^. 

Under the time-reversal symmetry, emergence of even 
number of Dirac points becomes trivial as discussed in 
Ref. 27: a Dirac point at fe = fep is always accompanied 
by its pair at fc = — fep- If there is only a single pair of 
Dirac points at fc = ifeg, they can be merged only at 
the special fc-points invariant under the transformation 
k -H> ~k, which are given as fc = G/2 with a reciprocal 
lattice vector, G. 



C. Generalized von-Neumann-Wigner theorem 

Most of recent studies on the point contacts mainly 
focus on ones discussed in the previous section, i.e., the 
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ordinary (two-fold) point contact in 2D lattices with two 
bands, which is located at a general fc-point. In this sec- 
tion, we develop a formalism for the general cases; the 
7n-fold contacts in d-dimensional systems with Ug bands 
(m < Ug), which are located at general or special k- 
points. We call it generalized von-Neumann- Wigner the- 
orem. 

Let us begin with the case of m ~ rig. The Hamiltonian 
at the Bloch wave vector, fc — (fci, fc2, ■ ■ ■ , fc^), can be 
represented as a m x m matrix. We expand it in the 
same form as Eq. (2): 



H{k) ^ Eo{k)i + R{k) ■ & 



(m) 



(11) 



where Eo{k) and the (m^ — l)-dimensional vector Jf2(fc) — 
(i?i(A;), i?2(fc), • • • , Rm^-i{k)) are real, and the traceless 
Hermite matrices, (t^™-* = ((Ti,(72,--- ,a„i^^i), denote 
the generalized Pauli matrices, which span the whole lin- 
ear space of the m x m Hermite matrices. For example, 
in case of m = 3, they can be chosen as the Gell-Mann 
matrices; 
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(12) 



If we find an m-fold contact of energy eo at fc = fco i there 
exists a unitary matrix U, which satisfies UH{ko)U^^ = 
eol- This implies H{ko) = eol- [m? — 1) conditions, 
R{k) = 0, should be satisfied at fc = fco, and the energy 
of the contact is given as eg ~ Eo{ko). 

Now, let us proceed to the case of rig > m. Instead of 
direct diagonalization of rig x rig Hamiltonian, H(k), one 
could renormalize it into smaller dimension following the 
formalism of Brillouin and Wigner^^. Let us divide the 
rig-dimensional Hilbert space S into m-dimensional sub- 
space Sa and (rig — TO)-dimensional subspace Sb- Then, 
the Hamiltonian matrix and its resolvent read 



^ ( HAAik) HABik) 

^ ' V^BA(fc) i?BB(fc) 



G(fc,e) 



H{k] 



( GAA{k) GAB(fe) 
^GBA(fc) GBB(fc) 



(13) 
(14) 



In particular, the resolvent matrix in the subspace Sa 
can be written as 

GAA{k,e)^(e-H^-''\k,e)y^ (15) 

with the energy dependent effective Hamiltonian 

H^''^\k,e) = Haa + Hab [e-HBBy' Hba- (16) 



The poles of GAA(fc, e) gives all the exact eigenenergies 
of H{k), as long as the projection onto Sa of their eigen- 
states do not vanish. 

As in Eq. (11), the effective Hamiltonian can be ex- 
panded as 

H(°ff)(fc, e) = Eoik, e)i + R{k, e) ■ d-^'^K (17) 

where Eq and R are real. Applying the previous argu- 
ment to this effective Hamiltonian, we can immediately 
see that a contact appears at fe = fco and e = eg, which 
is the solution of 



Eoik,e) 



and 



R{k,e) = 0. 

Equation (18) can be written explicitly as 

me - Tri?(<=^)(fc,e) 

= me - TtHaa - ^ Wb{e - (b)^^ = 0, 



(18) 



(19) 



(20) 



with Wb = {b\HBAHAB\b) > 0, where (b and \b) (6 ^ 
1, 2, • • • , rig — 2) denote the eigenenergies and eigenstates 
of -ffBB, respectively. Solving this equation with respect 
to e, we generally find (rig — m + 1) real solutions. 



which satisfy 

6 < Ci < 6 < C2 < 



+ 1), (21) 



(22) 



because the last term in the left-hand side of Eq. (20) 
diverges to ±00 at e — )- Cfi T and e — >■ ±00. If some of 
the eigenenergies of Hbb are degenerate, the solution ex- 
ists at that degenerate value, ^b+i = Cb = Cb+i- Strictly 
speaking, the number of solutions is reduced if there are 
some Wb = 0, namely the subspace Sa is disconnected 
from state \b). In such case, \b) is the exact eigenstate of 
the original Hamiltonian, H{k), and can be dealt sepa- 
rately. 

The m-fold contact takes place at fc = fco, which is a 
solution of 



Rj{k) = R{k,^j{k)) = 0, 



(23) 



and its energy is given as eq 
of equations is m^ — 1: i?,- 



^j(fco). Again, the number 



IjM — Rj,2 — ■ ■ ■ — ilj^m'-'^i 

First, let us consider the m-fold contact located at the 
general fc-points, which implies that the there are riu = 
d unknowns, i.e., (fci, • • • , kd). We further assume that 
there exist ric constraints near the fc-points where the 
contact occurs. In the vicinity of fc = fco, they should be 
expressed as 



0. 



s«.i?,(fc)=0 (r = l,2,... ,ne). 



(24) 
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within the hnear order of Rj{k), where Sj {i = 
1,2,- •• ,nc) is a hnear independent set of (m^ — 1)- 
dimensional vectors. The equation (23) is not ovcrdc- 
tcrmincd, if 

= "-u — fn^ + ^ + > 0, 

introduced in advance in Eq. (1), is fulfilled, and then the 
TO-fold contact of rid-dimension becomes feasible. Thus, 
we call Eq. (1) feasibility (generalized von-Neumann- 
Wigner) condition. 

Next, let us discuss the m-fold contact located at the 
special fc-points. In this case, the number of unknowns 
decreases: must be redefined from the spatial di- 
mension, d, to the manifold dimension of the special k- 
points. For example, 2D lattices with a refiection sym- 
metry have special fc-points invariant under the refiection 
(on the symmetry axis or on the boundary of the Bril- 
louin zone), which results in riu = 1- Instead, the number 
of constraints often increases at the special fc-points (See 
§ (11 F)). If the Eq. (1) holds after this redefinition of 
TT-d and TT-c, the equation (23) is not overdetermined, and 
the m-fold and rid-dimensional contact at the special k- 
points becomes feasible. 

We can imagine the cases that there exist the maxi- 
mum number of constraints (ric = (m^— 1)) on some spe- 
cial fc-points. Then, we always obtain i? = on those 
special fc-points, and the contact is classified into the 
essential degeneracy. Otherwise, the contacts occur at 
some unknown fc-points, and it is assigned as an acci- 
dental degeneracy. 



D. Dirac points in 2D multi-band systems 

This subsection is devoted to more detailed consider- 
ation on the ordinary point contact (rid = 0, m = 2) 
in 2D lattice with more than two bands. We assume 
that it is located at a general fc-point, fc = fcg (i.e., 
Uu = d = 2). Then, Rj.i{ki,k2) = -^^,2(^1,^2) ~ 
-Rj. 3(^11^2) — 0, is fulfilled at fc = fcp, and the num- 
ber of the constraints required for the feasible contact is 
one (ric = — + — 1 = 1). In the vicinity of 
k = ko, this constraint should be explicitly written as 

s,-R,{k)=0, (25) 

in the linear order of Rj{k), where Sj is a nonzero real 
3D vector. 

Now, the problem becomes equivalent to the Ug ~ 2 
one discussed in §11 B. Thus, we can derive some gen- 
eral results immediately. Dirac points, around which the 
energy bands always split linearly, are always stable, at 
which Xj = VkiRjlk^ko and Y, = Vfe^HjIfc^fco are lin- 
early independent. As derived in appendix. A, they arc 



evaluated within the linear order of Jfc = fc — fco as 

e±(fc) = £0 + D-^ (^A-Sk± ^{R, ■ Cy + DRfj (26) 
Ai = DVk.ij + • C, = DVk.^j + Yj ■ C, 
B = V^Eo, C = V,R, D = {I- Bf - >l. 

Here, all derivatives are evaluated at fc = fcg and e — eo, 
and Rj[k) is expanded as Rj = XjSki + YjSk2. The 
elliptic Dirac cone tilts if A = {Ai, A2) 7^ 0. We do not 
need explicit functional form of Cj{k) for the evaluation 
of Xj, Yj, and V^^j. In fact, they can be explicitly 
written as 

= Vfe,H+(l-B)-i(Vfe,So)C, 
Yj = Vk,R+ (1 - B)-' (Vk.Eo) C, 
VfcCj = (1 - By'VkEo. (27) 

Conversely, if the point contact at fc = fco is unstable, 
Xj and Yj are linearly dependent. In this case, there ex- 
ists Sq ^ satisfying Xjdqi + Yjdq2 = 0, and the band 
splitting at fc = fco -l- i5fc shows quadratic dependence 
on Sk parallel to Sq. The Dirac points are always cre- 
ated or destructed as a merged pair under the variation 
of the lattice parameters, which keeps the constraint of 
Eq. (25). Thus, also in the multi-band systems we gener- 
ally find even number of Dirac points between each pair 
of the adjacent energy bands. Under the time-reversal 
symmetry, the same argument with the one presented in 
the last paragraph of §11 B holds. 

In multi-band systems, one band can have Dirac points 
with both upper and lower bands, which we note upper- 
and lower-pair, respectively, since each Dirac point has 
its own pair. If the merging of upper-pair occurs at the 
same special fc-point with that of the lower-pair, it can 
be regarded as a new class of merging of Dirac points. 
Figure 1(b) shows the case where upper-pair starts to 
merge along Sq, while the middle and lower band starts 
to touch at the same time. The latter touching is induced 
by the merging of lower-pair along different Sq from the 
upper-pair one. We show in §111 A 2 the case where four 
Dirac points actually merge and form a three-fold point 
contact with one Dirac cone. 



E. Band overlap 

In §11 C we have shown that the Dirac points can be 
understood as stable feasible point contacts also in the 
multi-band systems. Now, let us assume that the Dirac 
points exist between the l-th and (/ + l)-th bands, and 
the system has a commensurate filling factor, v = 21. 
In order to see the interesting physics particular to the 
Dirac electron systems, the Fermi level should lie exactly 
at the Dirac points. Such a situation is realized, only 
when the l-th and {I + l)-th bands have no band overlap, 
i.e.; no hole or electron pockets. Otherwise, the Fermi 
level falls off the Dirac point, as shown in Fig. 1(c), and 
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the low energy excitation is dominated by the carriers in 
the pockets. 

At Tis ~ 2, the band overlap is absent when the di- 
agonal element, Eo{k), in Eq. (2) is the fc-independent 
constant, since e±(fc) = Eq ± |i2(fc)|. This condition 
is realized in a bipartite lattice, where all the diagonal 
elements (Eq and R^) arc fc-indepcndent. However, in 
7T.S = 3 the band overlap is possible even if all the diag- 
onal elements are fc-independent, e.g., when there is no 
direct hopping between sites with the same indices fi. 

For Jig > 3, we consider the two- fold contact (m = 
2) and again adopt the effective 2 x 2-Hamiltonian in 
Eq. (16). In analogy to Us = 2, it is straightforwardly 
concluded that the band overlap is absent if Eo{k,e) ~ 
iTri/^°^)(fc, e) is fc-independent (e-dcpendence is al- 
lowed), i.e., the solutions of Eq. (20), £,j (j = 
1, 2, • • • , rig — m + 1) (Eq. (21)), become fc-independent 
constants. This can be explained as follows; When the 
l-th energy band, ei{k), is equal to £,i, the l-th and (^ + 1)- 
th bands should have a contact at this fc-point. Thus, the 
upper bound of the l-th band is given by ^i. Similarly, 
we can also see that the lower bound of the (/ 4- l)-th 
band is given by These facts clearly show that the 
band overlap never occurs when ^^'s are constants. It 
should be noted that the fc-independence of £'o(fc,e) is 
a sufficient condition, and thus one has a chance to find 
systems without band overlap, even when Eo{k,e) is k- 
dependent. We show in the next section the examples 
of the geometry of lattices which could avoid band over- 
lap and could afford Dirac points at the Fermi level even 
without this sufficient condition. 



F. Symmetries and constraints 

In the present subsection, we consider the relation be- 
tween the symmetries and the number of constraints of 
Eq. (24). Since the spin-dependent cases are separately 
discussed in §11 G, we neglect the spin-degrees of freedom 
here. The following formulation is applicable to any type 
of symmetry present in the system. 

Let the Hamiltonian, H{k), be invariant under a sym- 
metry operation at a certain fc-point unchanged by this 
operation (invariant k-point). The effective Hamiltonian, 
Hi'^^)(k,e), keeps this symmetry, if Sa is chosen as an 
invariant subspace. Throughout this subsection, we con- 
sider only the 2x2 effective Hamiltonian (m = 2) unless 
otherwise noted. There are two possibilities, according 
to whether 7j(°^-'(fc,e) is invariant under the similarity 
transformation by a unitary operator, U, 



(28) 



or by an antiunitary operator, A, 

AH^''^\k, e)A-^ = i?(^ff)(fe, e). (29) 
Equation (28) implies the presence of a 2 x 2 unitary 



matrix, C/, which fulfills 

UH'^''^\k,e)U-^ = F(°*^)(fc,e) 



(30) 



Now, remind that any 2x2 unitary matrices can be ex- 
pressed in the form, 

U = e"^ (wo/ + • 6") , (31) 

where 0, ujq, uj ~ (oji, aj2, i^a) are real, and fulfill ujq + 
|a;p = 1 (since any unitary matrix can be written as a 
product of a phase factor and a SU(2) matrix). 

By inserting Eqs. (17) and (31) to Eq. (30), we obtain 
the commutation relation, 







(cff) 



(33) 



[u;-&,R-&]^2i{u;x R)- a, (32) 

which immediately gives 

u>xR = 0. 

Here, we used the formula, 

(A- a)(B ■ &) = {A- B)i + i{A-K B) ■ &. (34) 

Equation (33) indicates that u; should be either zero or 
parallel to R. When u; = (i.e., U = e"^l), it is trivial 
that the symmetry gives no constraint. As for u; 7^ 0, 
there should be two linearly independent vectors, s^^^ and 
perpendicular to a;, which serve as two constraints. 



•i? = 0. 



(35) 



at the invariant fe-points. 

Regarding A in Eq.(29), we always find a 2 x 2 unitary 
matrix, [/, which satisfies. 



UH''''^>{k, e)U-^ = ij('=ff)(fc, e). 



(36) 



since any antiunitary operators can be expressed as prod- 
ucts of complex conjugate and unitary operations. By 
using the expansion. 



(37) 



derived from Eq. (17), we can rewrite Eq. (36) as an 
anticommutation relation, 

|[/(T2, • crj = 0. (38) 

Then, by inserting Eq. (31), and using Eq. (34), we obtain 

0=1 (^loqI + iLj ■ &^ a2,R ■ o-j 

= |iti;2/ + s ■ R ■ o-| 



= 2iuj2R ■& + 2{s- R)I 
with s = (w3,wo,— wi), which gives 



UJ2R = 
s-R = Q. 



(39) 



(40) 
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Thus, we find a single constraint, s • il = 0, for a;2 = 0, 
and three constraints, R = 0, for u;2 7^ at the invariant 
fc-points. In general, the explicit matrix representation 
of U is required to judge which of these two cases is real- 
ized. However, if A is expressed as a product of a spatial 
symmetry operation, S, and time-reversal, T, one could 
know the number of constraints without it. Since S is 
commutable with T, we obtain = {STY, where the 
right hand side is represented by the matrix, UU* . On 
the other hand, 0^2 = is equivalent to U = , and 
thus UU* = i due to the unitarity, IJ-^ = = (^7^)*. 
These facts show that cj2 = is fulfilled when and only 
when = 1 holds. (Strictly speaking, = 1 must be 
fulfilled only within the invariant subspace at an invari- 
ant fc-point.) 

The representative example of a;2 = is the space- 
time inversion symmetry, i.e., the invariance under the 
spatial inversion, I, after the time-reversal, T, which 
satisfy 1^ = 1, identically. Both I and T give rise to 
the inversion of fc-points, k o — fc, and thus the general 
k-points are kept unchanged after the space-time inver- 
sion. Then, we obtain a;2 = 0, and a single constraint, 
s ■ R = at every fc-point. The 2D lattices with space- 
time inversion symmetry afford Dirac points at the gen- 
eral fc-points, because the feasibility condition for the 
point contact, = — + 1 + nc, holds there, for 
{nd,n^,m,nc) = (0,2,2,1). 

The example of ^2 7^ is provided by the invariance 
under the glide refiection, Q after the time-reversal, T. 
To be more concrete, consider the 2D lattice which is 
periodic in the x- and y-directions by lattice constants, 
a and &, respectively. If this lattice is invariant under 
the glide reflection in the x-axis, i.e., the translation by 
(a/2,0) after the refiection across the x-axis, the special 
fc-points on the Brillouin-zone boundary, = ±7r/a, 
are symmetry- invariant. In fact, the operation, QT, 
change that fc-points as (±7r/a, fcy) — > (^7r/a, — fcj,) = 
(in/a^—ky) — > (±7r/a, fcj,). At these special fc-points, 
— exp (2ikxa/2) = —1 gives 102 7^ 0. Actually, we 
obtain uj2 = 1, i-c, U = e"^(T2, because UU* = —1 indi- 
cates U = -([>-!)* = -(JJt)* = -U^. As a result, we 
obtain R = 0, and a line contact is found at the special 
fc-points, kx = 7r/a, which is an essential degeneracy. 

Sometimes, more than one invariances of Eqs. (28) 
and/or (29) are present at a certain special fc-point. Here, 
we consider the 7r-band of graphene, using a tight-binding 
model of the isotropic honeycomb lattice. Let us focus 
one of the unit cells, and assign the two carbon atoms 
in it as A and B-sites. At the K and K'-points, i.e., 
at the corners of the hexagonal Brillouin zone, Eq. (28) 
holds for both the reflection across the perpendicular bi- 
sector of the A-B bond, and the rotation by 2tt/3 around 
the A-site. By using the Bloch basis localized on A- 
and B-sites given in Ref. 12, the reflection is expressed 
by the orthogonal matrix, Ui ~ cti, since it exchanges 
A- and B-sites in the unit cell. By contrast, the ro- 
tation do not include the A-B exchange, and is repre- 
sented by the diagonal unitary matrix, U2 = diag(l, a) = 



Q!~"'^(cos(27r/3)i + i sin(27r/3)o-3), where a phase factor, 
a = e''^^^^, appears only for the B-site Bloch basis be- 
cause the B-site moves away to the other unit cell. These 
two symmetries requires lji x R = IjJ2 x R = with 
uji = (1, 0, 0) and 1^2 = (0, 0, 1), and thus R = 0. In this 
way, two point contacts at K- and K'-points in graphene 
can be understood as essential degeneracies. 

The above essential degeneracies can be understood as 
a limiting case of the accidental degeneracies. In fact, 
even when the geometrical anisotropy is introduced, i.e., 
in the absence of the reflection and the three-fold rota- 
tional symmetries"'^^, the honeycomb lattice continues to 
show point contacts at general fc-points which are the ac- 
cidental ones under the space-time inversion symmetry. 
Also, the decorated Honeycomb lattice we discuss shortly 
in §111 B shows the accidental point contact, in which the 
decoration site breaks both the space-time inversion and 
the rotational symmetries but keeps the reflection sym- 
metry instead. Therefore, the rotational symmetry in 
the isotropic honeycomb lattice only works to make the 
degeneracy essential. 

G. Spin-dependent cases 

Let us briefly discuss here the cases where the spin de- 
grees of freedom affects the orbital degrees of freedom, 
e.g., the systems under a commensurate spin-density 
wave (SDW) formation or with spin-orbit interaction. In 
the former example, we should consider the enlarged unit 
cell with SDW periodicity (magnetic unit cell) . The most 
general form of the Hamiltonian is 

where haa'(k,e) (cr, a' =t,i) are x ng matrices. 
The matrix dimension of the Hamiltonian is thus dou- 
bled. Here, refer to spin- dependent cases as those with 

%(fc)_7^0or/itt(fc)7^%W- . 

Special attention should be paid for the space-time in- 
version symmetry, because the energy bands show the 
Kramers degeneracy, which is the two-fold essential de- 
generacy seen at every fc-point (see appendix. C). We 
thus need to explore the four-fold degeneracy (to = 4) in 
order to consider the contact between two distinct dou- 
bly degenerate energy bands. In Refs. 34 and 35, the 
authors investigated the number of parameters required 
to describe the Hamiltonian with spin-orbit interactions 
in the context of topological insulators. Here, we reinter- 
pret their argument in our context. 

To discuss the four-fold contact, we choose two orbitals 
for each spin in such a way that 5'a is an invariant sub- 
space of space-time inversion. Adopting the method men- 
tioned in §11 C, we obtain a 4 x 4 effective Hamiltonian, 

[h[f^\k,e) h[f\k,e)) 
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The time-reversal operation, 7", changes the spin- 
dependent Hamiltonian as 



(a) bond-center (b) site-center 



(c) antisymmetric (d) mixed 



-(cff)* 

Hi 



(-fc,e) 
(-fc,e) 



while the spatial inversions as 



Uh)^^>{^k,e)U-^ 
Uh\f\-k,e)U-^ 



Uh^.f\-k,e)U-'' 
Uh^:P{-k,e)U-\ 



(44) 

with a unitary matrix U satisfying L02 = 0. Thus, the the 
space-time inversion symmetry, IT H^'^^\k^e)T^^I = 
H'-''^'>{k,e), implies 



uh 



(eff)< 

u 



{k,e)U- 



{k,e)^-Uh^'P^{k,e)U^ 



(45) 



Two spin-diagonal blocks, /ii'i"'(fc, e) 



and ft-M^' (A;, e), are 



/off) 

not independent, and we can expand the Hermite matrix, 
h'"^^\k, e), in the form of Eq. (2) with four real functions, 
i?o, Ri, R2 and i?3. Further, the spin-offdiagonal block 
is expressed in the form. 



4f (fc,6) 



{Zi{k,e)+tZ2{k,e))s-& (46) 



with two real functions, Zi and Z2, and the unit vector, 
s ~ {uj^^ujQj—LOi). Therefore, to find a contact, we only 
need to search the solution of Ri — R2 — R3 = Zi — 
Z2 ~ 0. This means that the number of conditions are 
reduced from — 1 = 15 to 5, and that the space- 
time inversion symmetry imposes 15 — 5 = 10 constraints 
on the Hamiltonian. To have feasible point contacts at 
the general fc-points, we still need extra three constraints 
in two dimension (riu = d = 2), because the number of 



1- 
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necessitated constraints is given as ric = ? 
by Eq. (1). 

As pointed out in Ref. 35, at the special fc-points, k = 
G/2, we sometimes have larger number of constraints, 
where i?('=*^)(fc, e) can become invariant under the time- 
reversal and inversion, separately. This case is discussed 
in appendix. B. 




FIG. 2: (Color online) Schematic illustration of (a) bond- 
centered and (b) site-centered inversions. Cross symbols indi- 
cate the inversion centers, circles and squares denote the lat- 
^ tice sites, and the shaded regions are the unit cells. Pairs of 
inequivalent sites which exchange by inversion are connected 
by broken lines. Pairs of inequivalent bonds connecting the in- 
equivalent sites which do not exchange by inversion are shown 
in solid lines. Panel (c) is the systems with asymmetric bonds 
whose Hamiltonian is invariant under the space-time inversion 
which transforms the phase of site-1 at the inversion center. 
Panel (d) is the mixture of bond-centered and site-centered 
sites. 



(a) ns=3 



(b) ns=3 



(c) ns=4 



(d) ns=A 



m — ■ — o 



(e) ns=4 






(f) ns=4 



(g) ns=4 




reflection axis 









— ti 


• 


1 ^ 


)Sa 


— t2 







— t3 








t4 


■ 
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3 ' 


)Sb 


= t3' 








= t4' 



III. TIGHT-BINDING MODELS 

In this section, we demonstrate that the feasibility 
(generalized von-Neumann-Wigncr) condition, Eq. (1), is 
useful to design Dirac points in unexplored lattices. We 
consider tight-binding model, which includes ng atomic 
orbitals per unit cell. Its Hamiltonian has transfer inte- 
grals between the sites and the potentials, W^, on the 
j/-th site (while we set W^, = unless otherwise noted). 
We basically consider single orbital per lattice site, but 
can deal with the multi-orbital ones, since the sites with 
I atomic orbitals can be regarded as I inequivalent sites 
at the same position. We introduce Bloch basis, fc). 



FIG. 3: (Color online) Representative lattice structure with 
bond-centered inversion of Us — 3, 4. Filled and open circles 
and squares represent the lattice sites, fi = 1,2,3,4, where 
1 and 2 belong to subspace Sa and others to Sb- Indices 
of bonds (ti) represented by solid and broken lines are given 
in the left panel. Panel (d) is the model lattice of a-ET2l3 
which is simplified to panel (c) by removing bonds with small 
transfer integrals. The lattice structure of panel (c) is reduced 
to that of panel (b) when site-4 is removed. If we take the 
vertical bonds in (f) uniform as, ts = t'^ = — t'4, we find 
the lattice in (g) which has glide reflection symmetry {yo/'2- 
translation in the vertical direction and reflection with respect 
to j/-axis) and affords essential line degeneracy. 
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m/ti = -1 



(tl, t2, f3)=(1, -0.5,1) 




kx>^ 



(9 ,<t))=(90, 180 ) 



(90, 90 ) 



FIG. 4: (Color online) Examination of the point contacts of 
the lattice structure in Fig. 3(b). (a) tijtx-tzlti diagram; 
the Dirac points are stable in the shaded region at W^j = 
0, (/J = 1 — 3). These stable regions shifts to those indicated by 
hatches when Ws/ti = 1 and -1, respectively, (b) Examples 
of band structures with Dirac points in the shaded region in 
(a), which are viewed from two different directions, {8,(j>) = 
(90°, 180°) and (90°, 90°). The polar angles, 9 and (j>, are 
defined in (fca;,fcj;,e)-space in the panel, where arrows indicate 
the direction of a view. 



(a) 



It3/tll 




(C) 

{tl, t2, t3, f4)=(-1 ,1,0,2) (-1,1,2,2) 



1 2 3 4 5 




(b) 



It3/tll' 



(1,1,0,2) 




2 3 4 
It4/tll 



FIG. 5: (Color online) Examination of the point contacts of 
the simplified a-type lattice structure in Fig. 3(d). |t3/ti|- 
|f4/ti| diagram for {a.)ti = —t2 and {h)t2=ti. Dirac points are 

a/(*4Ai 



stable in the shaded regions ((a 

\2 I /, /, \2 



ir-(i3Ai)2 < 2 

and (b): [ts/tif + {t^/tiY < 4). Panels in (c) show four 
different examples of the band structures with Dirac points. 



which has a probability amplitude only at the /z-th site in 
the unit cell i' — 1,2, ■ ■ ■ , rig). Then, the matrix ele- 
ments of Hamiltonian is given as Hi^i^{k) = k\H\v, k). 

The practical procedure to explore the fc-points is given 
as follovifs. Suppose we examine whether a given d- 
dimcnsional lattice structure could afford m-fold and rid- 
dimcnsional contact. Here, the lattice is characterized by 
a set of lattice parameter, t, which is dealt explicitly as 
variables. 

(i) Count the number of constraints, ric, on the general 



fe-point (riu = d). (If the special fc-points have 
larger number of constraints, separately count ric 
there, and redefine riu as the dimension of manifold 
formed by these special fc-points). 

(ii) Check whether the feasibility condition Eci.(l) is 
fulfilled. 

If fulfilled, we go onto the next step to specify the region 
of lattice parameter space with stable contacts, and to 
know the location of contacts in fc-space, fc = fcg. 

(iii) Start from the Ug x Hamiltonian, which is an 
explicit function of t besides k. Derive an effective 
m X m Hamiltonian, H^''^^k,e,t), using Eq.(16). 

(iv) Eqs.(18) and (19) give the analytical relationship 
between t and fco, when the contacts are present. 

One can thus identify the region with contacts within 
the lattice parameter space. Inside this identified region, 
the contacts are stable, whereas on its edge the contact 
becomes unstable. The examples given in this section 
follow these treatments. 

In the 2D lattices shown in Figs. 3, 6, and 8, the 
site indices /i = 1,2,3,4 are assigned to those repre- 
sented by filled circles, open circles, filled squares, and 
open squares, respectively. For the two-fold degeneracies 
(m = 2), we consider the 2D subspace, Sa, spanned by 
fi = 1,2. The m ^ 3 fold contact is also discussed for the 
Lieb lattice. §111 A and §111 B are devoted to the spin- 
independent systems, and the spin-dependent example is 
presented in §11 G. 



A. Space-time inversion 

In the present subsection, we specify the discussions on 
the space-time inversion symmetry to the tight-binding 
models, where the types of inversion symmetry are classi- 
fied by their appearance of lattices which directly reflect 
the shapes and locations of atoms in the crystal. We 
assume without the loss of generality that the atomic or- 
bitals have either odd or even parities under the inversion 
operation. Let us regard a pair of sites as equivalent or 
incquivalent, according to whether or not they are differ- 
ent only by a lattice vector. Inversion invariance of the 
lattice are basically classified into two cases: the bond- 
and site-centered inversions, which are shown schemati- 
cally in Figs. 2(a) and 2(b), respectively. If the inversion 
symmetry is accompanied by the exchange of the pair of 
incquivalent sites, it is called bond-centered one. In fact, 
the inversion center should locate at the middle point of 
the bond which connects that pair. On the other hand, 
the inversion symmetry without exchange of incquivalent 
sites is called site- centered one, since the inversion center 
can be placed at one of these sites. There is one more 
class of lattices which differ from the previous ones by its 
appearance, shown in Fig. 2(c), where the bonds exchange 
their sign (but not their absolute values) by the inversion 
operation. If the system includes pairs of incquivalent 
sites with and without the exchange (see e.g., Fig. 2(d)), 
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one can separately choose each pair as the Hilbert sub- 
space S'a, to which the bond-centered and site-centered 
inversion symmetry are adopted, respectively. Both cases 
will afford feasible contacts. 

In the bond-centered inversion, the Hilbert subspace 
S'a is spanned by |1,A;) and \2,k), where /i = 1 and 2 
are the indices for the pair of the corresponding atomic 
orbitals on the incquivalent sites exchanged by the in- 
version operation. In this case, the space-time inversion 
symmetry requires 

leading to {ujo,u)) = (0,1,0,0), where cti denotes the 
exchange between two incquivalent sites. Thus, wc im- 
mediately find a constraint, s ■ Rj ^ 0, with s = (0, 0, 1) 
at general fc-points. 

In the site-centered inversion, on the other hand, the 
subspace S'a is spanned by |1, fc) and |2, fc), where n — I 
and 2 are the indices of an arbitrary pair of incquivalent 
atomic orbitals which remain unchanged by inversion. In 
this case, we can see 

with {ojQjU}) = (1,0,0,0). Thus, wc have a constraint, 
s ■ Rj — 0, with s = (0, 1, 0) at general fe-points. 

The remaining issue is the inversion represented by 
Fig. 2(c). The subspace Sa is spanned by |1, fc) and |2, fe), 
where |l,fc) is connected to fc) (:^ = 2 ~ Ug) by the 
pairs of bonds which exchange their sign by inversion. 
Namely, the inversion operates in such a way that the 
phase of |1, fe) is shifted by tt as, 

^(cff) ^ 2:7-7j(off)7— ij-i ^ CT3i7(<=ff)*CT3. (49) 

Such case can be realized when |1, fe) and |2, fe) are spa- 
tially anisotropic and isotropic atomic orbitals (e.g. d- 
and s-orbitals), respectively. 

1. Bond- centered inversion 

Typical lattice structures with the bond-centered in- 
versions are given in Figs. 3(a)-3(f). All these lattices ful- 
fill the feasibility condition and could afford Dirac points 
in the certain range of lattice parameter space. 

Let us first examine the staggered square lattice in 
Fig. 3(b). We find that the Dirac points appear in the 
shaded region displayed in Fig. 4(a) on the plane of t2/ti 
and t^/ti at W^/ii = 0. The on-site potential, W3, which 
does not break the inversion, keeps the region with sta- 
ble Dirac points, as the figure shows for the cases with 
Ws/ti = ±1. The Dirac points in this region are indeed 
stable against the variation of lattice parameters includ- 
ing the on-site potentials. 

Next, we add one extra site to Fig. 3(b), which yields 
a lattice shown in Fig. 3(c). This lattice is a simpli- 
fication of a-ET2l3 shown in Fig. 3(d). The parame- 
ter regions which afford Dirac points in Fig. 3(c) are 



examined for ti = ~t2 and ti = t2 in Figs. 5(a) and 
5(b), respectively. The region with Dirac points of 
the former is extensive over the wide parameter space, 
\/ {ti/tiY — {t^/tiY — 2, which spreads along < 
1^3/^1! ~ 1^4/^1! < +00. Whereas, the region in the 
latter is confined to, (ta/ti)^ + (i4/ti)^ < 4, namely 
within the small value of lia/ti], 1^4/^1 1 < 2. It should be 
noted that parameter values of Q;-ET2l3 corresponds to 
the former case'^". Figure 5(c) shows the band structures 
for four different choices of lattice parameters. The last 
panel is the one on the edge of the Dirac point region 
in Fig. 5(b) (unstable Dirac point), where the two Dirac 
points merge and the dispersion in the merging direction 
become parabolic. 

Among the lattices with bond-centered inversions, the 
ones in Fig. 3(b), 3(c), and 3(f) have i?(<=ff)(fe,e) with 
fe-independent diagonal elements and thus have no band 
overlap. One can generalize these cases and refer to it 
as, "the effective bipartite lattices have no band overlap", 
where the words, "effective bipartite lattices" , mean that 
the lattices have no direct hopping between site-1 and 
-2. and at the same time, no indirect path, 1 — > — > 1 
or 2 — > ^ 2 (zy = 3, 4). In such case, the diagonal 
elements of Haa has no fe-dependence and the second 
term of Eq. (16) also has no diagonal matrix elements. 
Thus, the effective Hamiltonian indeed fulfills the suffi- 
cient condition to avoid band overlap. 

In the lower panel of Fig. 3(c), the sites = 3,4 and 
their related bonds 1^ u ^1 and 2— >■ i/ — )-2 are replaced 
by the the dotted lines, which schematically describes the 
effective bipartite lattice. This gives the physical inter- 
pretation of Eq. (16); Its second term formally describes 
the hopping of particle between site-1 and -2 in Sa me- 
diated by the occupation of particles at site-3 and -4 in 
Sb, and is regarded as having effective transfer integrals 
between 1— >-l and 2— >2. The resultant effective lattice, 
simplified from the a-type one, is topologically equiva- 
lent to the anisotropic square lattice in the lowest panel, 
which is known to afford Dirac points in a wide parameter 
region, as discussed in the context of a-ET2l3'^^. 

Finally, we briefly discuss the case with more than one 
invariances. Figure 3(f) has a bond-centered inversion 
and affords feasible contact points. If all the vertical 
bonds are taken uniform as = t'^ = t^ = t'^, one finds 
Fig. 3(g) which has the glide reflection symmetry (see 
§11 F for details) in addition to inversion. By the trans- 
lation of half the unit cell length in the vertical direction 
together with the reflection against the vertical axis, the 
site exchanges as 1 -O- 2 and 3 O 4. while the bonds 
remain unchanged. The special fc-points invariant under 
this glide reflection is fc^^ = tt. In fact, one can easily 
check that H^'^^^ spanned by Sa = {I72} automatically 
fulfills Eqs.(18) and (19) at = tt, which indicates the 
realization of essential degeneracy. The essential degen- 
eracy thus can often be understood as a limiting cases 
of feasible accidental degeneracy under the variation of 
lattice parameters due to the introduction of additional 
invariance, which we have also seen in the honeycomb 
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(a) ns=3 (b) ns=3 (c) ns=4 




FIG. 6: (Color online) Representative lattice structure with 
site-centered inversion of ns = 3, 4, (a) Kagome, (b) Lieb, and 
(c) anisotropic square lattices. The notation of lattice sites 
and the transfer integrals follow those given in Fig. 3. 



lattice of graphene in §11 F. 



2. Site-centered inversion 

The representative lattice structures with the site- 
centered inversion symmetry are the Kagome, Lieb, and 
the anisotropic square lattices shown in Figs. 6(a)-6(c). 
Noteworthy is the Dirac points in the Kagome lattice, 
which is found over the whole parameter region of the 
anisotropy of transfer integrals and on-site potentials. 
Under the variation of the transfer integrals, the Dirac 
points appear over the whole fc-space. Figure 7(a) shows 
the trajectories of Dirac fc-points between upper two 
bands for the several fixed values oi < t2/t3 < 2, under 
the variation of ti/ts. The values of ti/ta to have Dirac 
points along the M-F-X-M line are shown in Fig. 7(b). 
In the Kagome lattice, there always exist two sets of 
Dirac points between upper two and lower two bands 
at fc = ±fe(^ and ±fc(^, respectively. Here, these fc-points 
are described in the form, = {kiQi -\- ^292)72, and 
= [k^gi — fc2g2)/2, with the reciprocal lattice vectors, 
gi and 92 (sec Fig. 7(a)). Figures 7(c) and 7(d) show the 
energy bands for several choices of lattice parameters. In 
the regular Kagome lattice with isotropic transfer inte- 
grals (see the first panel of Fig. 7(c)), one pair of Dirac 
points between the upper two bands merge at F-point 
and touch the flat band in the parabolic manner. 

When one of the three bonds in the Kagome lattice is 
subtracted (e.g., is — > 0), the Lieb lattice is realized. In 
this limit, two Dirac points between the upper two bands 
and another two between the lower two bands merge at 
M-point (feo = {qi + 92)72), and a three- fold contact 
appears. Around the three-fold point contact, the center 
band is flat, and the upper- and lower- Dirac cones touch 
it at their tips. The merging of four Dirac points is shown 
in Fig. 7(e) in the vicinity of M-point. This particular 
class of band touching takes place regardless of the values 
oitiM^ 0). 

It is interesting to confirm that the feasibility con- 
dition, = — rn^ + ric = 0, in Eq. (1) holds for 
the three-fold point contact at general fc-point [n^ = 0, 
riu = 2 and to = 3) in the Lieb lattice. The Hamiltonian, 




FIG. 7: (Color online) Dirac points of the Kagome lattice, 
which merge at F and M-points, but do not open a gap 
throughout the whole (ti, t2, t3)-space, and under the on-site 
potentials, (a) Trajectories of point contacts between upper 
two bands in the anisotropic Kagome lattice in Fig. 6(a) un- 
der the variation of t^/ti for several choices of < t2/t\ < 2. 
(b) Variation of ts/ti along the M-F-K-M points which are 
indicated in bold lines in panel (a), (c) Variation of band 
structures along the F-M line with 12/ 1\ = 1. The last panel 
is the energy bands of the Lieb lattice, (d) Examples of band 
structure with and without on-site potential, (e) Merging of 
four Dirac points in the vicinity of M-point under t2/t\ — 1, 
ts/ti — > 0. Their cross sections along the different merging 
directions of upper- and lower-pairs correspond to those in 
Fig. 1(b). 



H{k), is a 3 X 3 Hermite matrix with only two nonzero 
elements, i/i2,-ffi3, which can be expanded as Eq. (11) 
with Gell-Mann matrices, Eq. (12). We immediately see 
i?3 = i?8 = from iJ^^ = 0, and Rq = = Q from 
H2J, = 0. Further, the space-time inversion symmetry 
offers R2 — R5 = R7 — 0. As a result, one finds Uc — 6 
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(a) 




(b) 



m/t 







\ A"honeycomb^ 





f'/f 




(e ,*)=(90, 60 ) 



FIG. 8: (Color online) (a) Decorated honeycomb lattice struc- 
ture of 71s = 3 without the inversion symmetry but with re- 
flection symmetry, (b) Upper shaded (Wa/t > {t' /tf /3 - 3) 
and lower hatched {Ws/t < —{t'/t)^ + 1) regions afford Dirac 
points between the upper two and lower two bands, respec- 
tively, and the region (shaded -I- hatched) at the center have 
the two sets of Dirac point, (c) Demonstration of the exis- 
tence of Dirac points under the reflection symmetry, which 
emerge along the K-K'-F lines in fc-space. The contacts on 
K- and K'-points at t' = O(isotropic honeycomb lattice) and 
those at t' /t 7^ are essential and accidental ones, respec- 
tively. The upper and lower pairs merge at F- and Af-points, 
respectively, at t' /t — ±3 and ±1. 



giving {ujo,uj) = (0,1,0,0), where the subspace Sa is 
spanned by |1, k) and |2, k). Thus, we can find two con- 
straints, s^^^-Rj = s(2).iJ^ 0, withs(i) = (0,1,0) and 
s(2) — (0^ 0, 1) on the special fc-points invariant under 
the reflection, located on the symmetry axis or on the 
boundary of Brillouin zone. 

We expect to find a feasible point contact on the spe- 
cial fe-points of symmetry axis, because the number of 
unknowns and constraints are = I, and 71^ — 2, re- 
spectively, giving nd = in Eq. (1). Figure 8(b) shows 
the parameter region with Dirac points on the plane of 
t'/t and Ws/t, whose origin correspond to the regular 
(isotropic) honeycomb lattice. The shaded and hatched 
regions afford Dirac points between the upper two and 
lower two bands, respectively. In Fig. 8(c), we show sev- 
eral examples of Dirac points at W3 = 0: When t'/t ~ 0, 
the essential Dirac points between the top and bottom 
bands are located at K- and K'-points at the corner of 
the hexagonal Brillouin zone. By the introduction of 
\t' /t\ > 0, these Dirac points split into two pairs and be- 
come accidental. With increasing |<'/t| they move along 
the symmetric axis connecting K-, K'-, and F-points. 
The lower pair merge first at M-point (unstable point 
at \t' /t\ = 1) and dissapear. Then, the upper pair merge 
at F point {\t'/t\ — 3) and dissapear. 



C. Dirac points in SDW 



in all, and sees that the feasibility condition is fulfilled: 
nd=:2 — 3'^ + l + 6 = 0. Although the Lieb lattice can 
have the three-fold point contact at a general fc-point, 
it always locates at M-point, a special fc-point, which is 
because it is generated by the merging of Dirac points. 

We finally mention that the anisotropic square lattice 
in Fig. 6(c) is effectively bipartite^ and thus has not band 
overlap. The Kagome and Lieb lattices also do not have 
band overlaps, even though the trace of iJ'^°*^)(fc, e) have 
fe-dependence and do not fulfill the sufficient condition 
to avoid band overlap. 



B. Reflection 

The final non-interacting example is devoted to the 
decorated honeycomb lattice shown in Fig. 8(a). Here, 
the inversion symmetry is broken by the decoration, 
namely an introduction of the third site described by 
square symbol, whereas, the reflection symmetry is re- 
tained, which exchanges /j, = 1, 2. 

In contrast to the space-time inversion, our consider- 
ation should be restricted to the special k-points on the 
symmetry axis, which is unchanged by the reflection. It 
requires 



We discuss the spin-density-wave (SDW) state as an 
example of the spin-dependent system, with LaOFeAs 
in mind. This material has five orbitals per site'^'^ 
and exhibits Dirac points in the vicinity of the Fermi 
level. Actually, they are reproduced within the self- 
consistent solution of the mean-field approximation of 
the five-orbital Hubbard model with on-site interactions, 
U, U' , and J, which denotes the intra-orbital-direct, 
inter-orbital-direct, inter-orbital-exchange Coulomb in- 
teractions, respectively^'^'^ . 

Here, we consider the simpler two-orbital Hubbard 
model on an isotropic square lattice, as shown in 
Fig. 9(a), and focus only on the (7r,0)-SDW, which has 
two- fold periodic spin modulation along the a;- direction^. 
The magnetic unit cell is doubled from the original one, 
and includes four orbitals. The spin-dependent mean- 
field Hamiltonian is an 8 x 8 matrix, and can be reduced 
to the 4x4 effective one described by Eq. (43). The 
space-time symmetry gives 10 constraints on the effec- 
tive Hamiltonian as discussed in §11 G. Further, there 
are two additional constraints, Z\ = Zi = 0, because 
the interaction processes do not mix states with oppo- 
site spins (j^^f^ = 0)- As a result, we can focus only on 



0-1 



the 2x2 matrix, h)^^ , to search for the contact, and 
the problem becomes equivalent to the spin-independent 
case. The only difference is that the space-time inversion 
symmetry is already used to reduce the problem, and 



(50) thus imposes no constraint on h. 



(off) 
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It is noteworthy that wc need more than one orbital per 
site in order to find Dirac points in the SDW state. To 
show this point, let us consider a single orbital Hubbard 
model (see Fig. 9(b)), where hff in the original mean-field 
Hamiltonian already has a 2 x 2 matrix form in Eq. (2). 
Since inter-orbital interactions are absent ([/' = ,/ = 0), 
the on-site interaction is simply evaluated as U n^^iif^]^ ^ 
U {{nf,f)n^^ + nf,^{nf,i) - {n^'f){Tifj,i)), where /i = 1,2 
denotes the index for the atomic orbitals in the magnetic 
unit cell, and a =t, i denotes the electron spin. This re- 
sults in i?3 = U {{nn) — (n24,)) /2 ^ in the expansion 
of h^f = EqI + R ■ a, and implies that the contact, i.e., 
the solution of i? = 0, is never found at any fe-point. 

Let us return to our main point, and consider the two- 
orbital Hubbard model, where the two atomic orbitals 
have dxz- and cJy^-symmetries as shown in Fig. 9(c). 
This model is invariant under the reflection across a;- 
or y-axis. We generate an effective 2x2 Hamiltonian, 
/i^'|.*^^(fc, e), focusing on two atomic orbitals (to = 2) in 
the original unit cell (^ = 1,2). Since the reflection ex- 
changes these atomic orbitals, 1 2, h^^'' = aih^^^^ai 
holds on the symmetry axis (rzu = 1), which leads to two 
constraints {ric = 2), i?2 = i?3 = 0. This situation is 
equivalent to the case discussed in §111 B, and thus the 
feasibility condition, Eq. (1), is fulfilled on the symmetry 
axis. Figure 9(d) shows one of the examples of self- 
consistent solutions of SDW by choosing the parameter 
close to those given in Ref. 33. The lowest two bands 
do show point contacts located on the symmetry axis 
{k, = 0). 



(a) 



magnetic unit cell 




(b) ns=2 



(ii,0)-SDW 



(c) ns=4 



^ (90, 270 ) 

hi 



dxz ([1=1) dxz (n=3) 2 
drz (|J,=2) dYZ (\i=4) 





FIG. 9: (Color online) (a) Square lattice with original unit 
cell and magnetic unit cell of SDW state. Arrows represent 
spins which have (tt, 0)-periodicity. Magnetic unit cell with 
(b) ?ia=2 and (c) Ws ~ 4. For the case with (c), Hamilto- 
nian is invariant under the reflection against either the xz- 
or j/x- plane, both of which exchange the two orbitals, dxz 
and dyz- (d) Energy band of the self-consistent SDW so- 
lution for parameters, t^fj, — — 0.2(nearest neighbor sites), 
(^'13,^24) ~ (0.3,0.15) (next-nearest neighbor sites in the 
(l,l)-direction), U = 1.2, U' = 0.9, and J = 0.15. The left 
panel is the cross-section at fc^: = as a function of ky, along 
which the Dirac point emerges (symmetry axis). 



IV. SUMMARY 

To summarize, we developed a general and simple for- 
malism to consider the feasibility of contact (degener- 
acy of energy bands), which can deal with the accidental 
as well as essential degeneracies. The cardinal stand- 
point of our framework is the feasibility (generalized von- 
Neumann-Wigner) theorem, = riu-^m^ + l — ric, which 
provides the number of constraints on the lattice, ric, 
necessary to have a feasible nd-dimensional contact with 
TO-fold degeneracy by some tuning of riu-unknown pa- 
rameters in multi-band systems. It enables us to judge 
without patiently solving the secular equation which lat- 
tice affords feasible contacts. Primarily, our framework 
provides a practical procedure to pick up only the degen- 
erate solutions of the secular equation, i.e., only fc-points 
with contact, selectively. This procedure plays an es- 
sential role in the design of Dirac systems, because the 
Dirac points often appear as an accidental degeneracy at 
unknown general fc-points. 

In simpler terms, the usual Dirac points at general 
fc-point in two-dimension correspond to the case with 
?id = (point contact), m = 2 (two- fold degeneracy), 
and Uu = 2 (fci and ^2), in which case, a single constraint, 
nc = 1 , is required to fulfill the feasibility condition. Such 



case can be practically explored in many actual crystals. 
In fact, besides a well known Q;-ET2l3, another 2D or- 
ganic crystal, (DIEDO)2-'4^ (X=Cl,Br), is found to have 
Dirac points whose space group is PI, i.e., only inver- 
sion symmetry present (for details of band structure, see 
Fig. 14 in Ref. 20). Even in the spin-dependent cases, 
the feasible contacts can be found in the same manner 
in the 2D systems, if some extra constraints besides the 
space-time inversion symmetry is present; we showed as 
an application to the LaOFeAs systems that the reflec- 
tion symmetry works as such extra constraint, and allows 
for the emergence of Dirac points in the spin-density- 
wave state. 

It is noteworthy that in three-dimension, the feasibil- 
ity condition for the point contact at general fe-point is 
fulfilled in the absence of constraints: Eq. (1) holds for 
(rid, ?iu, "1, ric) = (0,3,2,0). Therefore, the 2D lattice is 
favorable for the usual crystals which often have at least 
one symmetry such as inversion. 
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Appendix A: Derivation of Eq. (26) 

Expanding Eo{k^ e) — eo and i2(fc, e) within the hnear 
order of (5fc = fe — feg and 5e = e — eo, we obtain the 
self-consistent equation for 6e± = e±(fc) — eg 

5e± = (Sk ■ VfcSo + Se±S/,Eo) 

± \SkiVk^R + Sk2Vk,R+ Se±V,R\ , (Al) 

where all derivatives should be evaluated at fe = feg ^-nd 
e = eo- On the other hand, the self-consistent equation 
for 6£,j = £,j{k) — eo reads 



SS,j = Sk ■ \/kEo + S^j\7,Eo, 



which gives 



S^j = Sk ■ VfcO = (1 - ByHk ■ WkEo 



(A2) 



(A3) 



with B = VcE'o- The linearization of the relation, 
Rj{k) = R{k,£^j{k)), also gives 

Xj5h + YjSk2 = SkiVk^R+Sk2\7k2R + S^j^eR, (A4) 

where Xj = VkiRj and Yj = Vfe^iij. Thus, the self- 
consistent equation (Al) is rewritten as 

(1 - B)A± = ± \Xj6ki + Yj5k2 + CA±| , (A5) 

with A± = 6t± - 6ij = e±(fc) - ^j{k) and C V^R- Its 
solution reads 



A± = D-^ (Rj-C± J{Rj ■ Cy + DR 



(A6) 



where Rj{k) is expanded as Rj — XjSki + Yj6k2, and 
D = {1 — BY — is introduced. Then, we finally obtain 
Eq. (26), using e±(fc) = eo+(5^j + A± and = Sk-S/k^j. 
It is noteworthy that we do not need explicit functional 
form of £,j{k) to evaluate Xj, Yj, and Vfc^j. Actually, 
they are given as Eq. (27), owing to Eq. (A3) and (A4). 
Now, let us introduce the 2x2 matrix, 

M = {1 - B)i - C ■ 6- 

^ i + HABieo-HBBr^'HBA- (A7) 
Its eigenvalues are not less than one, since it satisfy 

mm = |||0)f + ||(eo-i?BB)"'^BA|0)f > 1, (A8) 
for any normalized 10) e ^a. Thus, we obtain 

{1- Bf -C^ = dot M > 1, (A9) 

and see that the value inside square root never become 
negative in Eq. (A6). Also, the factor (1 — B)~^ never 
diverges in Eq. (27), due to the inequality, (1 — B)^ > 1- 



Appendix B: Spin-dependent cases at k = G/2 

In Ref. 35, the spin-dependent cases are considered in 
detail at the special fc-points, k = G/2, where we can 
expect larger number of constraints than at the general 
fc-points. In fact, H^°^\k, e) is invariant under the time- 
reversal and inversion separately there, which leads to 



/ijf)(G/2,e) 



G/2,e) 



htr*iG/2,e) 

^(G/2,e), 



(Bl) 



and 



h'^:^,\G/2,e) = m:^,HG/2,e)U- 



(B2) 



respectively. In the case of a; = (i.e., U = e**^/), 
the number of constraints is unchanged. Otherwise, it 
is increased up to fourteen, since we obtain s^^^ ■ R = 
. i2 = and Zi = = from Eqs. (Bl) and (B2), 
where s^^^ and s^^^ are the linearly independent vectors 
perpendicular to w. In the latter case, a single extra 
constraint is necessary to find a feasible point contact on 
these special fc-points (n^ = 0), because Eq. (1) gives 
ric = rr? — 1 — riu = 15. However, such a feasible point 
contact is out of our main interest, because it is no longer 
accidental but essential. 

It should be noted that our interest differs from that of 
Ref. 35. Actually, the point contacts studied in Ref. 35 
is unfeasible in our context. It needs the fine-tuning of 
a single lattice parameter for its realization, or equiva- 
lently, an extra constraint for its feasibility. In terms of 
our formalism, they consider the point contact under the 



condition, 



1 



= -1. 



Appendix C: Kramers degeneracy 

In this appendix, we mention how the Kramers degen- 
eracy is derived from the space-time inversion symmetry 
using our formalism given in §11 C. Focusing on a sin- 
gle atomic orbitals, we obtain 2x2 effective Hamiltonian 
matrix, H^°^\ with the same form as Eq. (42), whereas 

^o-CT^'' are no longer matrices but scalars. We can adopt 
the consideration for the site-centered inversion, and ob- 
tain Eq. (45), which results in 



(Cl) 



because h^^E^ are scalars. Thus, the two-fold essential 
degeneracy takes place at every fc-point, since the effec- 
tive Hamiltonian always becomes a 2 x 2 scalar matrix. 
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